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An axisymmetric hydrodynamical model for the torus wind in 
AGN. II: X-ray excited funnel flow 

A. Dorodnitsyn 1,2 , T. Kallman 1 2 , and D. Proga 3 

ABSTRACT 

We have calculated a series of models of outflows from the obscuring torus 
in active galactic nuclei (AGN). Our modeling assumes that the inner face of 
a rotationally supported torus is illuminated and heated by the intense X-rays 
from the inner accretion disk and black hole. As a result of such heating a strong 
biconical outflow is observed in our simulations. We calculate 3-dimensional 
hydrodynamical models, assuming axial symmetry, and including the effects of 
X-ray heating, ionization, and radiation pressure. We discuss the behavior of a 
large family of these models, their velocity fields, mass fluxes and temperature, 
as functions of the torus properties and X-ray flux. Synthetic warm absorber 
spectra are calculated, assuming pure absorption, for sample models at various 
inclination angles and observing times. We show that these models have mass 
fluxes and flow speeds which are comparable to those which have been inferred 
from observations of Seyfert 1 warm absorbers, and that they can produce rich 
absorption line spectra. 

Subject headings: acceleration of particles - galaxies: active - hydrodynamics 

-methods: numerical - quasars: absorption lines - X-rays: galaxies 


1. Introduction 

One of the insights provided by observations of Seyfert galaxies and some quasars is 
the prevalence in their X-ray spectra of spectral lines and bound-free continua from ions 
of intermediate-Z elements . Early observations of Seyfert 1 galaxies using proportional 
counters and solid state detectors revealed spectra with strong absorption features in the 
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0.1-10 keV range (Halpern 1984). These features were attributed mostly to the edges of 
hydrogen and helium - like oxygen. The term ’’warm absorber” was proposed owing to 
the fact that the observed X-ray absorbing gas has an electron temperature lower than it 
would be if a similar level of ionization were produced by collisional ionization, ffowever, 
more detailed spectroscopic studies were hampered by the limited X-ray resolution of the 
ASCA and ROSAT satellites. The grating spectrographs on the X-ray telescopes Chandra 
and XMM-Newton provide unprecedented spectral resolution up to ~ 10 keV. These show 
that X-ray spectra obtained from ~ half of low-red-shift active galactic nuclei (AGN) contain 
many lines from ions of Fe, Si, S, O, Mg, and Ne , and that these are generally broadened and 
blueshifted by 100-500 km/s fjKaspi et al.l 120021 : ISteenbruggel 120051) . The presence of X-ray 
absorbing g as has been confirmed in the majority o f AGNs which are bright enough to allow 
detections f Reynold^ 1997 : McKernan et al. 2007 ) . The re is also a partial correspondence 
between UV and X-ray absorbers fl Crenshaw et al. 1999). 


X-ray observations of warm absorbers are consistent with the Seyfert 1/Seyfert 2 di- 
chotomy. For example, the properties of the X-ray emission in the Seyfert 2 galaxy NGC 1068 
co rrespo n ds to t he scatte red emission expected from warm absorbers in Seyfert 1 galaxies 
(Kinkhabwala et ah 2002 b 


Constraints on the position and dynamics of the X-ray absorbing gas can be deduced 
from the observed widths and virial argument s, and also from the variability studies of these 
spectra fjBehar et al.ll2003l:lNetzer et al.ll2003l) . These show an absence of correlated response 
of the warm absorber gas to rapid changes (~ days) of the continuum. This implies that the 
ionization time scale in the warm absorber gas is long (> months). Combined together, the 
line blueshifts, widths, and time variability analysis favors an origin of the warm absorber gas 
at R > 1 pc away from the Bid. This estimate coincides with the likely l ocati on of absorbing 
matter responsible for obscuration in Seyfert 2 galaxies flKrolik fe Begelmanl 11988 b The 


existe nce of an outflow from the torus has been s uggested bv iKrolik fe Begelmanl (11986 
1988 b and as the source of warm absorber flows by Krolik & Kriss ( 19951. 2001 b 


It is believed that this matter is in the form of a molecular torus which is responsible for 
obscuring the broad line region in Seyfert 2 galaxie_s, and which is thought to exist in most 
low and intermediate luminosity AGN (Antonucci & Miller 1986). A growing body of di- 
rect observational evidences advocates for the existence of the obscuring torus. Mid-infrared 
high spatial resolution studies of the nucleus of NGC 1068 using the Very Large Telescope 
Interferometer have resolved a dusty structure which is 2.1 pc thick and 3.4 pc in diameter 
flJaffe et al.l 2004 1). Observations support a multi-temperature model: the temperature of 
the warm component was established to be 300 K and inside of it a second, compact and hot 
(>800K) component has been found. Further studies of NGC 1068 systematically reduced es- 
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timates of the temperatures of different components (IPoncelet & Perrin 2006). Observations 
of the Circinus galaxy, which is among the closest prototype Seyfert 2 galaxies, also revealed 
a dense and warm T > 300 K componen t at about 0.2 pc from BH and cooler T < 300 K 
component at lpc fl Tristram et al. 2007). If the hotter component is located closer to the 
X-ray source, it may be attributed to the inner part of the torus, heated by the radiation 
of the compact nucleus. Although the evidence is strongest for nearby active galaxies, there 
is also a strong motivation to think that within the same obscuring torus paradigm exist 
those quasars whose central regions are heavily obscured by gas an d dust (Type II q uasar s). 
Evidence for this comes from spectro-polarimetric observations by Za ka m ska et ah (2006). 


This paper is part of a series whose main goal is to test the hypothesis that the torus is 
the origin for the warm abso rber flow. Preliminary results of this work have been reported 
in Dorodnitsvn et ah ( 2008 ) (Paper 1), in which we presented the results from a sample 
model and showed that the adopted model is promising in explaining the warm absorber 
phenomenon. In this paper we provide more details of our methods, and display results of 
models which span the space of input parameters. We present and discuss the hydrodynamic 
quantities which characterize our models: mass fluxes, velocity fields, and temperature struc- 
ture. We also show sample X-ray spectra, which we will discuss extensively in a later paper 
of this series. 


Our approach can be described as having three basic parts: i) setting up initial condi- 
tions, which requires defining an initial torus configuration and making assumptions about 
the external source of radiation; ii) implementation of the wind driving force (local heating- 
cooling rates and radiation pressure force) and actual 2D hydro dynamical calculations. The 
latter includes the numerical solution of the time-dependent 2D (so called 2.5D) system of 
equations, which takes into account centrifugal forces, and radiation pressure and heating 
terms; iii) calculating of the X-ray line spectra using a code which adopts Sobolev radiation 
transfer and ionization calculations for plasma in the intense X-ray held. Each of these steps 
is described in what follows. 


2. Governing equations 

We solve the following system of equations: 

( 1 ) 
( 2 ) 


f + V.(pv)=0, 


P 


<9v 

~dt 


+ (v • V)v = -Vp - pV<f> + pgrad, 
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— + V- (v(e + p + p$)) = H. (3) 

These are the conservation equations for: mass, momentum and energy. Heating and cooling 
processes are described by the function i7(ergcm _3 s _1 ); e - is the sum of the kinetic and 
internal energy densities: e = pv 2 / 2 + e. These equations should be supplemented by 
the equation of state which we assume to be polytropic: P = K p 1 , where 7 = 1 + 1/n, 
and n is the polytrope index and P = (7 — 1) e. A one-component, one-temperature T = 
Pp/ plZ, where p is the mean molecular weight per particle, 1Z = 8.31 • 10 7 ergK _1 g _1 is the 
universal gas constant and plasma with 7 = 5/3 is assumed to constitute the flow. All three 
components of the flow velocity v = (ty, vg, v</>) are calculated, assuming azimuthal ( d/dcf) = 
0) symmetry. Equations (HD-(J2D are cast in a non - dimensional form with the characteristic 
scales set by the properties of the plasma orbiting at a characteristic distance, Rq from 
a black hole of mass M§ ( in units of 1O 6 M 0 ). The characteristic scales are respectively: 
t 0 = R^ 2 f\jGM ~ 4.5 ■ 10 11 rpc 2 Mq 1 ^ 2 (s) for the time, where r pc is the distance in parsecs, 
and Vo = \/GM/R 0 ~ 6.6 • 10 6 M^ 2 r pc 1//2 (cms -1 ) for the velocity. 


3. Forces driving the flow 

Pleating and cooling of the gas 


The forces accelerating the wind in our model result from the gradient of gas pressure 
and from radiation pressure. The thermodynamic properties of X-ray heated gas depend on 
the spectrum of the incident radiation as well as on the local atomic physics. Under the 
assumption of photo-ionization equilibrium the thermodynamic state of photoionized gas 
can be parameterized in terms of the ratio of radiation energy density to baryon density 
( Tarter et al. 1969): 


where F x = L x e 

nr 

and r = 


£ = 47 t F x /n, (4) 

/ (47rr 2 ) is the local X-ray flux, L x is the X-ray luminosity of the nucleus, 

hip dr - is the optical depth, and n is the number density. We assume that the 
Jo 

attenuation is dominated by Thomson scattering n = 0.2(1 + Ah) — 0.4cm 2 g _1 , where 
A"h is the mass fraction of hydrogen, and the factor e -T , accounts approximately for the 
attenuation of the radiation flux on the way from the source toward a fiducial point. The 
methods adopted in thi s paper for treat ing the eff ects of radiation are essentially the same 
as those described in Proga et al (120001 ): IProgal (2007), and have been applied to various 
problems in the study of AGN and X-ray binaries. 
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Assuming that the there is a fraction / x of the total accretion luminosity Lbh available 
in X-rays and that the disk radiates a fraction T of its Eddington luminosity L edd = 1.25 • 


m 


10 44 M 6 we estimate: £ ~ 4 ■ 10 2 ■ f x T M e /(N 23 r pc ), where N 23 is the column density 
10 23 cm -2 . If the dynamical time within the flow is much larger than the characteristic 
time of the photoionization and recombination then the ionization balance is determined by 
the condition of photo-ionization equilibrium. The rates of Compton and photo-ionization 
heating and Compton, radiative recombination, bremsstrahlung and line cooling are then 
given by approximate formulas, modified from those of Bl ondinl (I1994T ) . for these processes: 


r IC (erg cm -3 s _1 ) = 8.9 ■ 1(T 36 £ (T x - 4T), (5) 

for the Compton heating - cooling; 

T x (erg cm -3 s _1 ) = 1.5 ■ 10~ 21 £ 3 / 4 T _1 / 2 (T X - T)T~\ (6) 

for the photo-ionization heating-recombination cooling , and for the bremsstrahlung and 
line cooling: 


A(erg cm -3 s _1 ) = 3.3 • 1(T 27 T 1/2 

+ (4.6 • 1CT 17 exp(— 1.3 • l0 5 /T)^ ( -°' 8 - a98 “ ) T-l/2 + 1(T 24 ) 5. (7) 

These formulae have been originally derived for a 10 keV bremsstrahlung spectrum (T x = 
2.6 • 10' K) and were found to be in a reasonable ( ~ 25%) agreement with numerical sim- 
ulations OBlondinl 19941 ) . Equations (}3|) - (ITI) are slightly modihed version of those of Blondin 
('1994b which accommodates new atomic data. Using the XSTAR code (IKallman & Bautista 
2001 ) we recalculated heating-cooling rates for the incident spectrum which is a power law 
with energy index a, and found results essentially equivalent to those given by equations 
0-0. Notice that in the case of a bremsstrahlung spectrum a formal value of a = 0 should 
be used in (JTj). For a power law with energy index a = 1.1 the results differ by < 30 % (see 
Figure [T]). Given these rates of energy deposition from the radiation to the flow, we write the 
total radiative heating-cooling function: H = Tic + T x — A. We have also performed several 
runs of our hydrodynamical models with different assumptions about heating-cooling, and 
found no important difference in the flow dynamics if using equations ([3])-(J7|) or the original 
formulae of Blondin (1l994h . and also between bremsstrahlung and power law spectra for sev- 
eral values of a. It appears that, for example, the effects of the optical depth are much more 
important. That is, the difference between curves for the power law and the bremsstrahlung 
spectrum at small £ (correspondingly high density) in Figure |T] becomes unimportant. 
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10 keV BREMSSTRAHLUNG a=1.1 POWER LAW 



LOG(T) 


Fig. 1. — Comparison of the results from the XSTAR X-ray photo-ionization code and 
analytic approximations o-©. Curves are plotted for different values of £ ranging from 
£ = 1 (lower curves) to £ = 10 4 (upper curves). Vertical axis: total radiative heating-cooling 
function H = Tic + T x — A in units of 10 -22 erg cm 3 s _1 . Horizontal axis: log(temperature). 
Curves; solid: XSTAR; dot-dashed: analytic. 
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The radiation pressure 

The radiation pressure force consists of the force due to continuum absorption g cont = 
F\j\k/c and due to lines: 


^rad = (■ Fjjyk/c ) M (t) , ( 8 ) 

where M(t) is the force multiplier ( Castor et al. 19751 ), and F uv is the local UV flux. We 
make use of the particular form (lOwock i et ah 1988): 


M{t) = kt a ((1 + r max ) (1 Q) - 1) /r^ a ", (9) 

where t = r/rj is the optical depth parameter, r] = Ki/a e is the line strength parameter, 
<7 e is the Thomson cro ss-se ction, and r ma v = trim ax . A parameter r/ ma x was introduced 
by Owocki et al. ( 1988 1 and Stevens fc Ka il ma n (1990) in order to limit the effect of very 
strong lines. That is, they assume a line number distribution which satisfies: dN/(drj dv) ~ 
7] a ~ 2 exp(— g/r/ max ), where N(rj, v) is the line number distribution. If / w x — > 00 , so that lines 


kt 


are distributed as a power law, one recovers the result of Castor et ah (119751 ): M(t ) 

In the opposite case of r max — > 0, the force multiplier is independent of t, and M max ~ fc?/ “ ax . 
As a result of this maximum line strength cutoff a correction factor appears in the relation 
for M(f), ([9]). The dependence of k and on £ has been numerically calculated and then 
fitted by the analytical formulae (IStevens fe Kallmanl 19901 ): 

k = 0.03 + 0.385 exp(-1.4£ a6 ), (10) 


r 6.9exp(0.16£°- 4 ), log 10 £<0.5, 

ogio^max | 9.1exp(-7.96-10- 3 £), log 10 £>0.5. 

From these, one can see that M(t) can depend sensitively on the ionization parameter. 
Taking a fiducial a = 0.5 (the value, adopted in all our calculations), one finds that M max = 
585 at £ = 0, then has two local maxima: M max = 724 at £ = 0.3, and M max = 743 at 
£ = 3.1. M max then drops to 1.7 at £ = 100 and decreases gradually to M max = 0.01 at 
£ = 1000 . 


3.1. Initial configuration: rotating torus with arbitrary Compton optical 

depth. 

We begin from a rotating toroidal configuration which is in equilibrium in the external 
gravitational field of the BH. The equation of state of the torus interior is described by the 


polytrope P = K p 1+1 / n . The distribution of the density (or pressure) in the torus interior 
was given by Papaloizou & Pringle (jl984j ) (in what follows refer to PP-torus for short), who 
assumed that the distribution of the specific angular momentum inside the torus is constant. 
In our case such a torus would not be in equilibrium because of the radiation pres sure from 
the central object. Thus we modify equilibrium equations oflPa palo i zou fe P ringle (I1984A to 
include the radiation pressure term. Since this cannot be done in a closed analytical form, 
we can write an approximate equation: 


p ~ i /i-r e - r w 1 

p n + 1 \ r 2 r 2 sin 2 6 

Note that equation (fTTl) must be understood as a bridging formula between two limiting 
cases: optically thin e -r ~ 1 (in which case it is the PP-torus with 1 — Y reduced gravity), 
and optically thick, when g ra d ~ 0 (PP-torus case). A constant C in the equation (fTIj) . 
parameterizes the distribution of the torus models and is connected with the distortion of 
the torus (this is described in more detail below). Including the radiation pressure reduces 
the effective gravity, and thus the torus gas needs less entropy to sustain it against vertical 
collapse. In both of these limiting cases this equation is exact. 



Notice that the problem of toroidal equilibrium in the presence of heating (or other 
radiation transfer effects) introduces a characteristic length scale through the optical depth 
r, leading to non-self-similarity of the model. Equation (ITTji was derived by assuming that 
the distribution of the specific angular momentum inside the torus is constant. Choosing 
non-dimensional units and working in terms of w, the cylindrical radius in units of R 0 , if 
we define the non-dimensional density p such that p{w = 1 ) = 1 , and the non-dimensional 
pressure P and internal energy e such that: P = (7 — l)e, and e(w = 1) = e 0 then 
n ( 1 -T) / 1 1 


e o = 


-C 


l+n 


The inner and outer edges of the torus are lo- 


n+1 (0.5 — C) n \x 2w 2 
cated at w~ and w + , respectively. Bounded configurations exist only for 0 < C < 0.5 and 
the distortion of the torus is described by the parameter d = ( tv~ + w + )/2 = 1/(2C). The 
boundary of the torus is matched to the exterior by the condition P = JT)~ 6 . The PP-torus 
is unstable to non-axisymmetric perturbations ( Papaloizou fe Pringle! 1984 1. However, this 
effect cannot be numerically investigated in the azimuthal symmetry which we adopt, since 
no signals can propagate in the (j) direction. At w > 1, matter that constitutes the torus 
has an excess of angular momentum with respect to the local ’equilibrium’ keplerian value, 
l(ru, z ) > Z(r), and vice versa, in the inner parts of the torus w < 1, Z(tX 7 , z) < l(r). It is the 
internal pressure of the torus (equation fill) that inhibits matter from settling to smaller (or 
larger, depending on angular momentum) orbits. The gas first evaporates from the part of 
the torus which is closer to the source of radiation and tends to settle at larger tv, as soon 
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as the back pressure supporting it drops. 

We begin our simulations from the stationary configuration determined from equation 
(ED- We follow the torus evolution as it is being heated by X-rays. No replenishing of 
the gas which constitutes the initial torus is provided: Therefore, the torus will eventually 
lose all its mass and will completely evaporate. However, in the regime we are looking 
for, the evaporation is not dramatic and does not significantly deplete the torus during the 
characteristic dynamical time. 


In the following sections we will show that the existence and character of the flow from 
the heated torus depends critically on the geometry. That is, it depends on the divergence of 
the flow streamlines, the strength and incident angle of the X-ray illumination, and on the 
direction of the effective gravity in the rotating frame of the torus. The flow is intrinsically 
two-dimensional, and therefore cannot be ad equately described a priori by ID models, such 
as those preformed by Chel o uche fe Netzer (12005 1. Furthermore, the shape of the torus, 
and thus the launching surface for the flow, is affected by_the f low. So the torus interior 
cannot be considered as a boundary condition (e.g. as in 
to include it in the computational domain. 


Balsara fe Krolik (1992J)) ; we need 


4. Methods 


For our computations we adopt a spherical-polar coordinate system ( r,9 ), extending 
the computational domain {ri,9j} from r m = 0.01 to r out = 50 in radius, and from 0 to 
7 r in the polar domain making no assumption about equatorial symmetry. The number of 
points in the radial, N r , and polar, Ng, directions are taken to be equal: N = 140, in low 
resolution and N = 300 in high resolution grids. The {r*} grid is non-uniformly spaced, i.e. 
r 2 = n + (r ou t - r- in )(kl /( - Nr ~ 1) - 1 )/((kr r/{Nr ~ 1) - 1)), and r i+1 = r t + (r* - n-i)k 1 r /{Nr ~ 1 \ 
for i — 2,N r » 1, and the refinement factor is k r = 4. In order to achieve better resolution 
of the flow itself rather than the torus interior we also adopt a polar grid which has non- 
uniform spacing S9 t = 9i — 0j_i, so that the maximum refinement is approached at 9 = 7t/4: 
S i+ i = Si/kl^ Ne 1 ' ) at 0 < 9 < 7t/4, and 9 = 7t/4: <5 i+1 = dikl^ Ne ^ at 7t/4 < 9 < 7t/2 (and 
analogously spaced in the southern hemisphere). Boundary conditions are axially symmetric 
at 9 = 0, 7T and outflowing at r in and r out 


ve numerically th e sys tem of hydrodynamical equations (00)- (E]) we use the code 
Stone & Norma nlll992 ). Note that the characteristic time of X-ray heating/cooling 


To so' 

ZEUS2D l 

can be much shorter than the dynamical time, which in such a case introduces strong stiff- 
ness to the system of equations (CO)-©- To overcome this difficulty, some modifications have 
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been made to the code. The most important one is the implementation of a fnlly implicit 
update of the energy in equation (j3J) just prior to the transport step in ZEUS2D. Addionally, 
we account for the radiation pressure (equation (J8J)) term. As an initial test we have evolved 
a toroidal distribution of matter for two rotational periods and found the configuration to 
be stable. The gas is illuminated by the incident X-ray radiation with a power law spectrum 
with an energy index, a = 1. The heating/cooling rates are described by the approximate 
analytical formulae give in equations 0-0. 


Warm absorbers 


We test output of our hydrodynamical models against the ability to predict warm ab- 
sorber spectra. To do this we use the output from the hydrodynamical code, of p, v, and T as 
an input to the calculation of X-ray line and photoelectric absorption spectra. The numeri- 
cal code has been specifically developed for calculation of spectra in the X-ray domain and 
makes use of procedures developed for the XSTAR flKallman k, Bautistal 120 01) code, while 
calculating the ionization structure and distribution of opacities, and treating the radiation 


transfer in the Sobolev approximation (Rybicki & Hummer 1983). Although the goal of this 


paper is to show that pure hydrodynamic 2D models can produce warm absorber spectra, we 
present here only sample spectra, assuming pure absorption. We postpone a more detailed 
discussion, including a full 3D transfer calculations, to a separate publication. 


Results. 


The most important parameters which determine the properties of the warm absorber 
flow are the initial Compton optical depth ry = t(9 = 90°) of the torus (or equivalently the 
maximum initial torus density n max ), and the distance from the BH, R 0 . We also explore 
the dependence on T and d. Other parameters are chosen having some typical values: the 
mass of the black hole: Mbh = 10 6 M 0 , the Compton temperature of the X-ray radiation 
T x = lOkeV, and the fraction of X-rays and UV radiation f x = fuv = 0.5. (For rotating 
flows exposed to a multi-temperature radiation, see e.g. Proga et ah ( 2008h ). The latter 


is consistent with typical energy distributions of the radiation close to the BH (ILaor et al. 


19971 ). We neglect any changes in the BH luminosity. The important thermal time scales 
within the flow, namely the Compton heating and cooling time, f x and the dynamical time, 


tdyn, may be of the same order t x ~ fy 


yn 


10 10 s. This is discussed in more detail later 
in this section. Thus the outflowing gas may not be in thermal equilibrium and adiabatic 
losses are likely to be important. Notice that a nearly hydrostatic Compton heated corona 
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can exist only at r < Ric = GMbu /i/7 ZT X ~ 8- 10 16 (M 6 /T Xj 7) cm, where T X)7 is the Compton 
temperature in terms of 10 'K. In all of our models the major flow is located at r » Riq. 

We have calculated 20 models, including combinations for: = 1.3 (models Aj), and 40 

(models Bi); R 0 = 0.5, 1, 1.5; T = 0.1, 0.3, 0.5 all with d = 2.5; and two models with d — 5 
(models C(). These are summarized in Table 1 where some of the characteristic results from 
the computed models are presented. In what follows we describe in detail the cases which 
best illustrate the most important results. We also discuss the dependence of our results on 
parameters, based on the behavior of the ensemble of models. 

The model A§ is similar to that described in Paper 1, although the initial torus in the 
model which is described here has a different distribution of p and e (see equation (JTTj) , and 
the discussion thereafter), and smaller t^. In Paper 1 this model has been described in 
detail. Calculations presented here reveal more details and confirm the conclusions of Paper 
1. We begin here by describing results from model Bq, and later discuss how it differs from 
model A (i . 

Model, B 6 has r/ = 40, R 0 = 1 and V = 0.5 and corresponds to a Compton thick 
(r/ ~ 40) torus having large n max = 10 7 cnf 3 and mass M tor = 9.3 ■ 10 5 M 0 . Results 
are displayed in Figure [2l where the evolution of the distribution of density is shown as a 
function of time (the density scale is such that 0 corresponds to 10' cm -3 ); Figure [3] where 
the distribution of pressure is shown at t — 3 (the pressure scale is such that 0 corresponds to 
4.7 • 10 -4 dyn/cm 2 ) ; Figure [I] for various quantities as a function of the inclination d; Figure 

[5] (left panel), where the effect of the distortion parameter d is demonstrated; and Figure 

[6] showing horizontal ’slices’ of the velocity and temperature at constant height, z. In the 
case of this model, the torus column is high enough to effectively screen the torus interior 
from penetrating X-rays. This leads to formation of a nearly pure funnel flow, i.e. the torus 
interior, and hence the shape of the surface responsible for launching and collimating the 
flow, is essentially unaffected by X-ray heating on time scales < f rot . 

Here and in what follows we discuss the time evolution of our models in terms of t, 
measured in units of the characteristic time of rotation, t 0 - After t = la high pressure 
region created by X-ray heating extends to r ~ 4.5 — 5 pc throughout the area that is not 
shadowed by the high density torus. At this time the torus is located at 6 < 50°. The 
distortion parameter has a value d ~ 2.5, i.e. the torus shape is almost unchanged from its 
initial value. This is shown in the upper left panel of Figure[2l Within the part of the flow 
which is not shadowed by the torus, high temperature gas expands in a spherical bubble with 
radius, r < 5.2 pc in which the temperature is T ~ 3 — 10 T vir (r), where T vir = 2.6- 10 5 M 6 /r pc 
is the local virial temperature. An axisymmetric region exists between w < 0.75 pc and z < 2 
pc where the temperature, T ~ 10 T vir (r) . That is, high temperature, T ~ 3 ■ 10 6 K, but 
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low density gas fills the torus funnel. The ionization parameter (equation [4|) in this region 
is £ ~ 10 4 — 10 5 . The outer edge of the torus extends to ~ 4.25 pc in temperature, and to 
~ 4.5 pc in density contours. 

Figure [2] (lower left) shows density and velocity fields for model Bq at t — 3. Figure 
[3] shows that a high pressure region expands to height z ~ 6 pc from the equatorial plane. 
The torus inner edge is inferred from the temperature and density maps to be w~ ~ 0.83 pc. 
Inside the of this radius, which we refer to as the torus throat, the temperature is T ~ 
10 6 — 10' K. A wide nozzle with (w max — w min )/z max — 2.12, where z max ~ 0.4 pc is formed, 
having inner radius of w ~ 0.85 pc. The torus outer edge is slightly shifted to w + ~ 4.5pc. 
The values of (the minimum ionization parameter along a radial line) and the column 
density vary significantly with the inclination angle. Figure 0] shows the distribution of 
radial and poloidal velocity, £, density and the rate of growth of number density with radius 
as function of 9 at t — 3 for model Bq. Near the axis, £ m j n (0 — 20°) = 10 4 and the 
column density is N23 = 10~ 3 . Note that if t x > tdyn, he. the gas is not in thermal 
equilibrium, then £ is not as meaningful as when f x << tdyn- When t x > td yn adiabatic losses 
strongly affect the temperature of the gas. At larger 9, the ionization parameter decreases: 
£ min (# ~ 25°) = 3-10 3 and at higher inclination, £ gradually reduces from £ min (# = 45°) = 12, 
eventually becoming £ mi n {6 ~ 60°) = 2.5. At a critical angle, 6 ~ 40°, a strong rise of the 
column density reflects the fact that the line of site penetrates the dense torus body rather 
that through the wind (c.f. Figure [4l lower right). The column density increases, from 
A23 = 0.3 at 9 = 45° to N23 ~ 100 at 9 ~ 60°, providing total obscuration. Figure [3] 
also shows the position of the sonic surface determined by the relation v p /c s = 1, where 
v p — ( v r + v e ) 1 ^ 2 is the poloidal velocity and c s = (7 ZT / fi) 1 / 2 is the speed of sound. Behind 
the torus a low entropy region exists which is bounded from the sides by a quasi-stationary 
shock. The existence of this structure can be understood from the following considerations. 
If the flow were perfectly symmetric in both hemispheres, then it should have v z = 0 at 
z — 0, and the z = 0 plane would be the equivalent of a rigid wall (reflecting boundary). 
Thus, if v z < 0 behind the torus the formation of a shock structure is anticipated. Generally, 
this is the kind of picture one expects to observe from a supersonic wind flowing over a rigid 
obstacle. 

At t — 5 in model Bq (Figure [21 lower right), the density maximum is located at 
w ~ 2 pc. The inner edge of the torus does not shift significantly from the position it has at 
t — 3: w~ ~ 0.75 pc in density maps (and ~ 1 pc in temperature maps); the outer edge is 
at w + ~ 4.3 pc. The temperature of the torus interior is in the range 10 3 — 6 • 10 4 K. A hot 
flow is located near the axis, bounded from the sides by the torus throat, and having high 
temperature: ~ few ■ 10 6 K. A significant drop of ionization parameter £ from ~ 6 ■ 10 3 to 
~ 6, occurs again at 9 > 45° — 50° (c.f. Figure SJ) , where the column density also rises from 
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IV 23 = 0.04, to N 2 3 = 30 at 9 > 60°. The asp ect ratio of the torus is: A = Rn/H ~ 1 i n 
accord with what is inferred from observations (IKrolik fc Begehn aril 1 19861: IJaffe et al.ll2004j) . 
At low inclinations, 9 < 10°, everywhere in the wind the poloidal component of the velocity 


is determined by v r . However, at 9 > 50° inside the torus throat, the vg component is 
important, i.e. vg ~ v r at w < lpc. 


Model Aq has Tp = 1.3, Rq — 1, T = 0.5 and is very similar to the model described 
in Paper 1. It differs from model B e in that the smaller optical depth of the torus interior 
cannot shield the gas from a extensive X-ray heating and the torus loses mass from large 
parts of its surface. The initial maximum density of the torus is n max = 10 6 cm -3 corresponds 
to initial torus mass, M tor = 9 • 10 4 M 0 . Figured shows the distributions of poloidal velocity, 
£, density and the rate of growth of number density with radius as a function of 9 at t = 3 for 
model Aq (in the same format as Figure [4]). During the evolution, a region of high pressure 
extends from r ~ 4.5 — 5 pc at t = 1 to r ~ 12 pc at t = 3, and to r ~ 20 pc at t = 5. The 
inner edge of the nozzle shifts slightly from w~ ~ 0.8 pc at t=l to w~ ~ 0.83 pc at t=3, and 
w~ ~ 0.75 pc at t=5. At later times the behavior of the model Aq is similar to models A\ 
and , and can be inferred from Figured 


It has been mentioned that in model B$ much of the torus interior is opaque to pene- 
trating X-rays. Remarkably, the minimum nozzle cross-section doesn’t change much at late 
times, implying that the mass-loss rate becomes quasi-saturated. Note that in the case of a 
ID flow M is roughly set by the position of the sonic point, which in turn is set by gravity. 
In the case of a 2D nozzle, the mass-loss rate is determined by X-ray heating, gravity and the 
minimum nozzle cross-section. In the case of model B$ the latter remains almost unchanged 
in time. We believe this model is probably most representative in showing the key features 
of X-ray excited flow. However, models A* may generally have broader angular patterns in 
which a warm absorber spectrum is observed, as will be discussed below. Only comparing 
synthetic spectra with observations can answer the question of what model is more adequate 
in describing the phenomenon of warm absorbers. 
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Model 

rl 

Rq 

T 

d 

mo° 

c 'max,t=3 

,,45° 

( -'max,t=3 

?; 10 ° 

max,t=5 

,,45° 

max,t=5 

M t 

=3 

M 

t=5 

Ai 

1.3 

0.5 

0.1 

2.5 

516 

155 

624 

332 

4.09- 

1-^ 

O 

1 

6.54 

1-^ 

O 

1 

CO 

A 2 

1.3 

0.5 

0.3 

2.5 

710 

317 

847 

330 

1.48 ■ 

O 

1 

CO 

4.31 

0 

1 

CO 

A 3 

1.3 

0.5 

0.5 

2.5 

707 

267 

760 

291 

2.34- 

0 

1 

CO 

2.14 

O 

1 

to 

A 4 

1.3 

1 

0.1 

2.5 

547 

189 

514 

217 

1.76- 

0 

1 

CO 

1.66 

0 

1 

to 

A 5 

1.3 

1 

0.3 

2.5 

526 

179 

605 

343 

9.65- 

0 

1 

CO 

6.34 

• 10" 2 

A6 

1.3 

1 

0.5 

2.5 

570 

235 

670 

337 

2.02 • 

10" 2 

1.8 

10" 2 

a 7 

1.3 

1.5 

0.1 

2.5 

360 

197 

413 

230 

6.21 • 

h-L 

0 

1 

CO 

1.64 

• 10” 2 

As 

1.3 

1.5 

0.3 

2.5 

388 

169 

540 

310 

1.40 • 

10" 2 

5.68 

• 10" 2 

Ag 

1.3 

1.5 

0.5 

2.5 

317 

207 

663 

370 

2.66- 

h-L 

0 

1 

to 

1.23 

0 

1 

Bi 

40 

0.5 

0.1 

2.5 

673 

318 

522 

320 

3.56- 

0 

1 

CO 

7.38 

0 

1 

CO 

b 2 

40 

0.5 

0.3 

2.5 

590 

257 

1004 

471 

1.16 ■ 

0 

1 

CO 

1.87 

• 10” 2 

b 3 

40 

0.5 

0.5 

2.5 

907 

383 

957 

459 

5.49- 

0 

1 

CO 

2.56 

0 

1 

to 

b 4 

40 

1 

0.1 

2.5 

506 

205 

438 

236 

3.04- 

0 

1 

CO 

1.53 

0 

1 

to 

b 5 

40 

1 

0.3 

2.5 

536 

216 

587 

276 

8.71 • 

h-L 

0 

1 

CO 

2.85 

• 10” 2 

B 6 

40 

1 

0.5 

2.5 

641 

271 

676 

324 

1.55- 

10" 2 

7.24 

• 10~ 2 

b 7 

40 

1.5 

0.1 

2.5 

395 

179 

496 

187 

1.39- 

h-L 

0 

1 

to 

3 • 

0 

1 

to 

b 8 

40 

1.5 

0.3 

2.5 

541 

185 

610 

329 

2.22 • 

10" 2 

7.17 

• 10” 2 

b 9 

40 

1.5 

0.5 

2.5 

547 

248 

602 

347 

3.2- 

10“ 2 

8.53 

• 10” 2 

Ci 

40 

0.5 

0.5 

5 

890 

464 

770 

349 

3.46 ■ 

0 

1 

CO 

1.16 

• 10” 2 

c 2 

40 

1 

0.5 

5 

789 

788 

772 

770 

1.35 ■ 

10" 2 

8.01 

0 

1 

CO 


Table 1. Models, for different initial T£, Rq, T, and d and results for the maximum velocity, 
^max T=time(^ ms_1 )) where $ is the inclination angle; and the mass-loss rate, MT=time(M 0 /yr). 


Mass loss within the funnel flow 

It is instructive to consider the distribution of variables within a horizontal cross-section 
at a certain height above the equatorial plane. In so doing, we interpolate the solution 
from an (r,0) - spherical grid to a (z,t u; 100x100) Cartesian grid. Figure [6] shows the 
distribution of temperature and z-component of velocity, in terms of the escape velocity, 
U e sc = (2GM B h/t) 1//2 , at different heights for model f? 6 . 

A hot region extends to tu ~ 1 pc at z = 0.2, and to w ~ 2 pc at z = 1. The ’’funnel” 
can be seen in distributions of both temperature and velocity. At the X-ray heated boundary 
of this nozzle gas is being heated so that its temperature increases suddenly to 10 6 - 10 7 
K. This fact reveals an analogy between the torus flow with X-ray excited winds in X-ray 
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binaries flBasko et al.lll977l: iMcCrav fc Hatchettl Il975h : we discuss this further later in this 
section. Notice that in our case the inner surface of the torus both serves as a copious source 
of a gas and as a collimating funnel. 


Figure [5] shows models B 6 and C 2 at t = 4 and Figure [8] shows density and velocity 
streamlines for models A\ and A 5 at t — 4. Notice, there is little difference between Figure 
|8] (left panel, Model Ai) and Figure [8] (right panel, Model A 5 ); the effect of smaller R 0 is 
partially compensated by the fact that T is also smaller, thus reducing the effective gravity. 
If T ~ 0, then uj~ ~ 0.5 (for C = 0.2 in the equation (1111) ). However, when V = 0.5, as 
in model Bq (Figure 0 left panel), the effective gravity at the innermost optically thin edge 
of the torus is reduced by half. Figure [5] (right panel) shows a model with initially large 
distortion d = 5 (C = 0.1), model C 2 in Table 1. 


In model A 6 , a well-developed wind is observed in the vicinity of the high density torus, 
following the equal pressure contours; the maximum radial velocity is observed close to the 
axis at v max (9 ~ 3°) = 700 km s -1 . As a general trend at t — 3 the maximum velocity has 
a plateau at 20° < 9 < 50°, v max = 220 km s -1 , and lower values closer to the equatorial 
plane (Figure [7]). The flow is approximately symmetric in both hemispheres. At later times, 
t = 4 and t = 5, the behavior of the model is similar to t — 3: namely, v max (9 ~ 4°, T = 
5) = 900kms _1 , and on the plateau being u max ~ 380km s -1 . The torus is losing mass in 
all directions, although with very different speed at different inclinations. Because we are 
solving equations of ideal hydrodynamics (with only a small numerical viscosity), accretion 
through the inner boundary (at r=0) is negligible: M in (M 0 yr _1 ) < 10~ 8 . The maximum 
mass flux per unit solid angle yr -1 sterrad^ 1 ) peaks at 9 ~ 13° at t — 3, i.e. at 

much higher inclinations than v max , and at 9 ~ 55°, M^ iax = 0.01 at t=l, = 2 • 10 -3 at 

t=3, and M^ ax = 0.02 at t = 5. The total mass-loss rate at t = 3 is M(M 0 yr^ 1 ) ~ 7 • 10~ 3 . 

The mass-loss rate is M(M Q yr _1 ) ~ 2.4 • 10~ 2 , at t = 4, and M(M 0 yr _1 ) ~ 4 • 10~ 2 , 
at t=5, and the change of the mass-loss rate with time is dM/dt((M Q yr _1 )/yr) ~ 10~ 6 . 
Comparing distributions of v and n we conclude, for example, that the apparent minima 
of v p ~ v r correlate (with a certain lag) with maxima of n and vise a versa, reflecting 
conservation of mass flux. 


As in model Aq, the model Bq funnel wind carries mass flux which doesn’t change much 
during the evolution. The maximum velocity is as high as ~ 1000 kms -1 near the axis, and 
typically 200 < w max < 600kms _1 at 15° < 9 < 50°. The bulk of the gas, which potentially 
may produce warm absorber features, moves with comparable speed. However the largest 
observed velocity in B 6 model is v max (9 ~ 3°) = 1200kms _1 , at t=5. The mass-loss rate is 
M(M 0 yr _1 ) ~ 3.4 • 10~ 3 , at t=3, and M(M 0 yr -1 ) ~ 7 • 10~ 2 , at t=5. 
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Spectra 


Computing absorption spectra is a key test for the warm absorber flow model. Several 
sample spectra are shown here, although the detailed discussion of methods and results of 
calculations of such spectra is postponed to a later paper. 

Figure [9] shows the model A§ spectrum observed at different inclinations. This figure 
shows the warm absorber spectrum at t — 3 and at t — 4. At t — 4 a rich X-ray line 
absorption spectrum exists in the range 43° < 9 < 52°, and in the range 47° < 9 < 55° at 
later times, t = 5. 

At t — 3 the Bq model predicts a rich spectrum for 42° < 6 < 47°. At later times 
a similar spectrum appears at lower inclinations. Figure flOl shows the model Bq spectrum 
observed at different inclinations at t = 4. At t — 5 the spectrum exists between 45° < 
6 < 50°. Notice that the region of the funnel wind in this model is bounded by the area 
unshadowed by the torus: 0° < 9 < 40°. At 0 > 30° column density becomes N 23 — 0.45 and 
the ionization parameter is £ < 20. At higher inclinations the X-ray flux in the 1 < E < 2 
keV range becomes severely absorbed. 

Figure [IT] shows the evolution of the observed properties of the warm absorber flow with 
time (in the same time units) for model A 3 . It can be seen that warm absorber spectra are 
changing slowly on a timescale At ~ 1. This is typical for most of our models and shows 
the range of times over which our solution can be considered as a representation of a steady 
state warm absorber flow. 


A quantitative analysis of our synthetic spectra and comparison with observations will 
be done in a later paper. This is due in part to the need for full 3-dimensional treatment of 
the transfer and scattering of line photons, which we do not present here. Rather, the spectra 
in figures [9], HQ] and HU are calculated assuming pure absorption. We can calculate crudely 
some of the properties of individual lines, and show that these are generally consistent 
with observations. A convenient way to do this is to discuss the profile of what is likely 
to be the strongest line in any synthetic spectrum, the La line of OVIII. In model A§ at 
t = 3 the full width at half-maximum (FWHM) of this line is ~ 200km s _1 . Closer to 
the BH, the maximum observed velocity is greater, i.e. models A 3 and B : > give FWHM 
~ 400km s^ 1 at 6 ~ 43°, and 40°, respectively. The centroid energy of the line is at a 
blueshifted velocity (50 — 200)kms _1 with respect to line center. These velocities are less 
than th ose observed from, eg., NGC 3783, but are comparable to those observed from other 
objects (IMcKernan et al.ll2007l) . Such comparisons should also include the effects of scattered 
emission, which may skew the line centroid and red edge, and which we have not considered 
here. 
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Analytical estimates of the mass-loss rate 


The mass-loss rate found from numerical calculations is in approximate agreement with 
theoretical expectations. The value of the mass-loss rate, M, can be estimated by integrating 
the average mass flux (j) over the surface area of the torus exposed to X-ray radiation, 
X ~ 27T 2 R 2 /A. where A = R n / H ~ 1. (j) may be est i mated using the same arguments as 
those of iBasko et al.l (119771 1 and iMcCrav fe Hatchettl (119751) . Namely, heating from a BH 
creates a narrow transition layer, a ’’skin” on the surface of the torus. There, temperature 
rises almost discontinuously from inner ’’cold” ( T ~ 10 4 , T < T vir ), to outer ’’hot” (T > T vir ) 
value. This transition can be seen in Figure [HI 

Matching momentum, p + pv 2 and mass flux, j = pv below and above this discontinuity, 
we obtain a well known relation: j 2 = (. Ph — Pf)/(pf l — p/ 1 ), where subscripts 0 and h 
refer to values below and above the discontinuity. Being heated, the gas expands and its 
specific volume, V = 1/p increases. Above the discontinuity the flow is assumed to be 
isothermal so that P ~ 1/V. In the P-V plane, the transition between points Pq, Vq and 
Ph, Vh goes through the straight line with an inclination, (P h — P 0 ) / (14 — Vo) > ( dP/dT ) T , 
and it follows j 2 < — {dP/dV)T = p\(? sh , where c s = (7 ZT/ p) 1 ^ 2 is the velocity of sound. 
Since j = PhVh , it follows ^ vy, < c s .h and the flow im mediately above the discontinuity is 
subsonic ( Basko et al. 1977 : McCray fe Hatchett 19751) . From the momentum conservation, 
P 0 ~ P m = P h + phi’l, and the mass flux associated with such heating, can be estimated as 
p m 

(j) = — -A — — — where M.h is the Mach number above the discontinuity, P m is the 

v h Mh{l + Mp 

pressure below the discontinuity, and for simplicity we assumed Vh — c s ,h- 

McCray & Hatchett (1975) have calculated the state of the gas in the optically thin 
layer of a stellar atmosphere heated by X-rays. From their results it follows that the relation 
between P m and F x can be cast in the form: P m = 10~ 12 a_i2 F x , where a i? ~ 1, reflecting 


the shape and effective temperature of the incident spectrum ( IBasko et al. 1977). Although 
it is essential (in order to obtain stationary transonic flow, correctly matching boundary 
conditions at infinity) that the flow above the discontinuity is subsonic, we assume that the 
sonic surface is located not far from the discontinuity, estimating i >2 = c s> h, M-h = 1- Next, 

we write: F x = - (1 + A ), where A is the effective X-ray albedo of the X-ray heated skin 

and we take A = 0.4 (which we simplistically assumed to be optically thin), and assume 
/i = 0.5. Calculating M =< j > X, we finally obtain: 

/x r m 6 


M (M 0 /yr) ~ 0.16 




A 


( 12 ) 


where is the temperature above the discontinuity in units of 10 6 K. Inserting relevant 
parameters, such as V = 0.5, / x = 0.5,i? = lpc, M 6 = 1, A = 1, and adopting the value 
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of Th taken from our numerical model Aq, Th — 10 6 K, we estimate the mass-loss rate: 
M(M 0 yr _1 ) ~ 4 • ICC 2 . Comparing results from this approximate formulae with those 
summarized in the Table 1 we conclude that they are in good accord. Given the torus mass 
in the A e model M tor = 9.3 ■ 1O 4 M 0 , we conclude that it may sustain such mass loss for 
~ 1 • 10 6 yr. The upper limit may be inferred from Table 1, and is found to be ~ 10 8 yr. 


Adiabatic loses 

The characteristic time scale at which the energy is deposited to the flow via Compton 
processes, t x can be cast in the form: 

Ms) ^ 9.4 ■ 10 1 " (13) 

where T = T/T x and T x = 2.9 ■ 10 7 K. This should be compared with the dynamical time, 
tdyn °f the flow: 


tdyn(s) ~ 4.3 • 10 i0 r pc Vf. 


(14) 


When tdyn tx, the outflowing gas departs from thermal equilibrium and one must account 
for adiabatic losses, A^, when calculating the temperature of the gas. Notice that the 
properties of the two-phase (or multiple-ph ase) gas are conventionally described by the S- 
curve on the T — E diagram (Krolik ct al.l ll98lh . where 5 = F x /(nkTc) is the other form 
of the ionization parameter. That is on the T — S plot those places where dT/dE > 0 
are stable to isobaric perturbations. Places where dT/dE < 0 are unstable. Including A a( j, 
may significantly lower the temperature of the hot phase f Chelouche fe Netzer 20051 ). This 


a e n 


T x to 


temperature can be estimated by equating the Compton heating rate, A I0 — 4/cF x 

m e c z 

V 

the adiabatic losses rate, A a d ~ —pc 2 . The flow near the funnel walls is less divergent than 

r 

it would be in the case of a spherically-symmetric wind, in which case the latter expression 
is a factor of 2 larger. Assuming that above the discontinuity v ~ c Sl we obtain: 


T h ( K) ~ 5.7 • 10 6 


f x T\ 2/3 


L pc 




(15) 


which gives Th ~ 2 • 10 6 K, for parameters adopted in this paper. This value is in good 
agreement with the value of Th, which is found from T{w,z) distributions shown in Figure 
E Three major regions within the funnel flow may be emphasized: i) a ” discontinuity” 
where temperature is rising from the inner ’’torus” value to Th ~ 10 6 K; ii) a ’’plateau” where 
T ~ Th and thermodynamic characteristics of the flow result from the interplay between Aj 0 
and A a( q iii) region of hot, over ionized flow where T — > T x . 
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Returning current 

From Figures ([5]) and (jBJ) we see that there exists a region, behind the dense torus, 
where outflow is switched to inflow. This gas rejoins the torus in the shadowed region. For 
example, taking the model A 6 , and integrating the mass flux over the region where v p < 0, 
we obtain M in = 4 • 10 -6 M 0 yr -1 at t=3. That makes ~ 6% of the total accretion rate 
M accr , required to maintain 0.5L ec jd luminosity of the BH, given the efficiency of accretion, 
rj = 0.06. At the same time, much more mass, ~ 2 • lO _3 M 0 yr _1 , is lost within the funnel 
(6 < 50°) in the X-ray excited wind. Matter that is removed from the funnel is replaced 
by gas from the torus interior. Thus, a weak large scale convection flow is observed in the 
simulations. This effect is most clearly seen in models with large d, such as model C 2 , shown 
in Figure [5] and is due to a strong drop of v p as the outflowing gas is passing the shock wave 
front behind the torus (c.f. Figure [3]) and being unable to escape from the potential well. 


Radiation force 


The dependence of the radiation pressure on the ionization parameter, £, is determined 
by equations (jHJ) and (HUD . In the region of the fast flow the wind is too overionized for 
the radiation force to be important. This resembles Low Mass X-ray Binaries (LMXB) 
case in which the radiation pressure is also found to be insufficient to drive a significant 
outflow (Proga & Kail maul 2002 ). The ionization parameter drops below ~ 100 at 6 > 40° 
- the value determined by the torus aspect ratio, A. Thus the radiation pressure may be 
of importance at higher 6 and at these inclinations its relative strength is determined by 
the attenuation of the X-ray and UV fluxes. For model A§ we have r ~ 6 at 6 ~ 90° and 
r ~ 1 at 6 ~ 60°; i.e. the torus becomes Compton thin at 6 < 60°. The radiation pressure 
exhibits complicated behavior with varying 0, having multiple maxima and minima. The 
force multiplier, M(t,£(6)) peaks at 44° at r ~ 2, where (j'rad/s'grav ~ 5. Generally, two 
maxima of g ra d are observed at a given 6 along a radial line. The second peak becomes smaller 
at higher inclinations, i.e. in models A; the radiation pressure is determined mainly by the 
properties of the X-ray heating (i.e. £(#)) rather than by the attenuation of the UV flux. At 
higher 9 smaller maxima occur at smaller r; the inner skin of the torus exerts considerable 
radiation pressure, although at large 9 it is opposed by the back pressure of the torus interior. 
We calculated a model which has the same parameter values as model A 6 but with g rac 1 = 0. 
At t — 3 this model gives v max (9 = 10°) = 564 km s -1 and v max (0 = 45°) = 194 km s -1 . 
Comparing with Table 1 values we see that for the range of angles where warm absorber 
flow is observed the radiation pressure doesn’t play a major role in the flow acceleration. 

In models ip the attenuation is much stronger than in models A* and consequently the 
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secondary maxima of g ra d which were observed in models A* are suppressed by the e~ T 
attenuation. The radiation pressure is important only on the skin of the torus but almost 
everywhere points in the wrong direction, opposing the back pressure of the torus interior. 
Only at 6 ~ 45° it points in the direction tangential with the torus surface, but as r ~ 3pc 
the density drops and £ rises so that M(t,£) becomes small. 


Dependence on T, -R c ,0;Ti> and d 


If the interior of the torus is optically thick to X-rays then the torus loses mass m ostly 


from the surface, much as in the ’self-excited wind’ scenario for X-ray binaries (I B asko et ah 


1977). As shown above, in such a case the torus throat serves as a funnel and the gas is 


injected to the flow from the funnel walls. 


Notice that the location of the narrowest part of this funnel determines the characteristic 
terminal speed of the wind. In order to explore this, we have made a set of runs similar to 
models A 3 and Aq, but with reduced For model A 3 which has T£ = 2 and Rq = 0.5, we 
find that for 6 = 10° and t = 4.5, the maximum velocity v™ ax equals 738 km s _1 . For model Aq 
which has rf = 2 and Ro = 1, we find v™ ax = 432 km s' 1 for the same 6. If in the latter model 
we make optical depth smaller, = 1 we obtain: Vp liax (T = 4.5, 6 = 10°) = 400 kms -1 . 
This shows, in accord with our expectations of the mass flux conservation, that the torus is 
losing mass from deeper inside. As shown in Table 1, reducing Rq has the effect of increasing 
the maximum velocity. An increase of T has the same effect. However, this maximum 
velocity may be observed at a different inclination. Increasing the distortion parameter, d 
has an effect of some increase of the maximum velocity, redistributing v™ ax {6) to higher 
inclinations. From numerical solution we notice that the torus aspect ratio, A = r/H ~ 1 
does not strongly influence the evolution. That is because it is the most inner part of the 
throat which determines the dynamics of an evaporative flow. This inner throat is located 
at high 6 so that it remains optically thick most of the time. Numerical experiments confirm 
that the geometry of this inner throat remains approximately unchanged in time. 


Figure [2] shows that the geometry of the innermost part of the torus, i.e. the densest 
part (roughly located between 0.5 and 2 pc) shrinks considerably in the vertical direction 
during the process of the evolution. This is the result of the joint action of the radiation 
pressure and the back pressure of the hot evaporative flow. This is particularly interesting 
as it resembles the geometric ally thick outskirts of AGN ac c retion disks which are know n 
to be unstable to self-gravity (IKolvkhalov fe Svunvaevlll980l: IShlosman &; Begelmanlll989l) . 
The physics of such systems is complicated, and is subject to various possible competing 
effects. The self-gravitating instability may operate also in the torus body, perhaps leading 
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to a d ynamical system of molecular-dusty self-gravitating clouds (as in Krolik & Begelman 
(119881) ) . If this is the case, the optical depth of the torus, is crucial as in the optically 
thin case th e torus will effectively cool and collapse to a thin disk with subsequent star 
formation ( Toom re 1964). In the other extreme (rf >> 1) the released energy can go 
to in crease the velocity dispersion of the clouds, effectively supporting the torus thickness 
(IPaczvnski 1 9781 ). Strong IR radiation pressure exerted on these clouds, which can come 


from internal reprocessing of X-rays, can produce significant vertical force ([Thompson et al. 


20051 : iHonig fc Beckertl 120071 ). and may suppress the sel f-gravity inst ability and at the same 


time provide pressure support against vertical collapse (IKroliki 2007). Vertical support, and 
partial suppression of gravitational collapse, may also be provided by radiation pressure 
from star formation within the torus or the obscuring flow (JWada & Norman 2002). Further 
heating and loosing mass induces a torus to expand and change of shape. 


We have calculated models Bq and A$ with 100x100 resolution further in time to learn 
the late time evolution. At T = 17 in the B 6 model the torus has two extended lobes in 
both hemispheres with an opening angle 45°. They have a certain degree of asymmetry with 
respect to the equatorial plane. The shape of the obscuring structure no longer resembles the 
initial torus; the column densities are in the range of N 23 = 10 at 30°, and N 2 3 = 10 3 at 30°. 
The radial velocities in this structure are in the range of 200 — 400kms _1 . Our model does 
not allow for the replenishing of the torus; obviously the torus will evaporate completely if 
given enough time. Thus, in the the model Aq the torus evaporates completely by the time 
T = 15. These results imply that in order to get a quasi-stationary warm absorber flow the 
replenishing time should be of the order of the mass-loss time. The whole torus configuration 
may be unstable in a secular sense; the instability is driven by the long characteristic time 
of the global torus heating/cooling (due to expansion, winds, radiation loses), advection of 
heat in the torus body by internal flows etc. For example, the mass-loss rate, M ~ £; the 
surface area E increases during the torus expansion. If after some time of extensive heating, 
the torus separates into several parts, further mass loss will increase due to the larger total 
surface area of the fragments. 


6. Conclusions 

We have studied X-ray excited winds from the putative gas-dusty torus in AGN. We ap- 
proach this problem using numerical methods combining detailed hydrodynamical modeling 
with calculation of the warm absorber spectra. Hydrodynamical calculations include two- 
dimensional, axially-symmetric rotating flow, driven primarily by X-ray heating. Compton, 
bremsstrahlung, and photoionization heating/cooling processes were taken into account as 
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well as the radiation pressure force, which was calculated in the Sobolev approximation. A 
code combining XSTAR for photoionization calculations with the Sobolev radiation transfer 
has been developed for the calculation of the spectra. 

We find that a rotationally supported torus heated by radiation from the inner accretion 
disk and black hole can indeed be a source of the material we observe in the warm absorber 
flow. We find that the inner throat of the torus is not only important as a source of the gas 
but also because it creates a funnel for the outflowing wind. This leads generally to larger 
velocities within the funnel, and different velocity distribution within the warm absorber flow 
from those derived from models based on spherically-symmetric winds. The wind mass-loss 
rate within the funnel is not very sensitive to the details of the initial torus distribution and 
approaches ~ 0.02 — 0.09 M 0 yr -1 . Strong X-rays heat the gas within the funnel, producing 
a fast, ~ 1000 kms -1 , ionized flow near the axis, and slower, < 500 km s -1 , flow closer to 
the funnel walls. This is where optical depth effects become important and a warm absorber 
spectrum is produced. Using methods developed in studies of X-ray binaries we were able 
to estimate the mass-loss rate from such funnel flow, finding it to be in a good agreement 
with our numerical solution. 

The funnel flow is found to be promising with respect to obtaining high velocity warm 
absorber flows. What is beyond the scope of our models is the possibility of having mul- 
tiple phases in such high velocity flow, on spatial scales smaller than our grid resolution. 
Our treatment of the gas thermal properties will produce two-phase behavior at our grid 
resolution; we do not find this behavior, owing to the fact that the cooling timescales are 
generally too long. The answer to the question of whether there can be high velocity ’bullets’ 
or ’embedded clouds ’ on length scales smaller than the resolution of the grid is related to 
the problem of the origin of broad and narrow UV/optical line emitting clouds, and requires 
different computational methods from those employed here. 

Our models which have initial Compton depths r 0 > 1, aspect ratio Ro/H ~ 1, and 
located at 0.5 < r < 1.5 pc predict warm absorber spectra, thus confirming the main 
conclusion made in paper 1. The existence of such spectra depends on the fact that the 
flow is intrinsically two-dimensional, meaning that both the dynamics of the funnel flow is 
different from ID models and optical depth effects are important as they strongly depend on 
inclination. The latter point requires that we include the entire torus in the computational 
domain rather than considering it as a boundary condition. The distribution of the ionization 
parameter, £ depends strongly on 6 , further confining the range of angles where conditions 
are right for the warm absorber flow to be observed. In most of our models warm-absorber- 
like spectra are produced in a 10° range, at 6 ~ 40 ± 5°. This range is set both by the 
initial aspect ratio of the torus, which we take to be ~ 1, and by the thickness of the X-ray 
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heated ’skin’ of the torus. More optically thin models produce warm absorber-like spectra for 
9 ~ 40 ± 10° , as they potentially provide more partially optically thin gas for evaporation. 

The bulk of the gas in this scenario has a terminal velocity of the order of the escape 
velocity at the inner torus edge. Because of the funnel mechanism part of the gas is re- 
distributed to lower inclinations and acquires a higher terminal speed, ~ 1000 km s 1 . In a 
real AGN environment such flow may contain clumps and irregularities and even dust, which 
are not captured in our studies because of the intrinsic limitations our methods. Accounting 
for the multiple phases of a gas (on a subccllular level) may reveal this in more detail and 
may also broaden the range of angles where the warm absorbers appear. 


The part of the flow that is shielded by the optically thick part of the torus body can 
also flow out as part of a torus global expansion. Thus it strongly depends on the deposition 
of energy directly to its interior. This problem is related to oneyrf the infrared support of the 
AGN torus vertical structure against gravitational collapse ( Kr olik 20071) and also requires 
additional investigation. 
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Fig. 2. — Model Bq] Color- intensity plots of the logarithm of the dimensionless density as 
a function of time which is given in orbital periods. In the northern hemisphere this is 
superimposed with velocity vectors. Axes: distance in parsecs. 
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LOG PRESSURE AT t=3 
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Fig. 3. — Model i? 6 ; Color- intensity plots of the logarithm of the dimensionless pressure at 
t = 3; at northern hemisphere superimposed with velocity vectors. The location of the sonic 
surface is marked in red. Axes: distance in parsecs. 




Fig. 4. — Model Bq ; distributions of v r , and the poloidal velocity v p (upper left), £ (upper 
right), n (lower left), and the rate of growth of number density with radius (lower right) 
at time, t = 3. Curves are marked by a inclination angle 9. Horizontal axis: distance in 
parsecs. 
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Fig. 5. — The effect of distortion parameter d; Velocity streamlines superimposed on contours 
of the number density. Model Bq (left panel), model C '2 (right panel) at time, t=4. Axes: 
distance in parsecs. 



Fig. 6. — Model Bq at time, t = 3. Horizontal ’slices’ of temperature and z-component of 
velocity are marked at each curve by the value 2 in parsecs. Vertical axis; left: z-component 
of the velocity in terms of local escape velocity; right: temperature. Horizontal axis: distance 
in parsecs. Curves; dashed: velocity; solid: temperature. 
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Fig. 7. — Model A 6 ; As in Figure [4] here are shown the distributions of v r , v p (upper left), 
£ (upper right), n (lower left), and the rate of growth of number density with radius (lower 
right) at t = 3. Curves are market atop by an inclination angle 6. Horizontal axes: distance 
in parsecs. 
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Fig. 8. — The effect of T; velocity streamlines superimposed on contours of the number 
density. Model A\ (left panel), model A$ (right panel) at time, t=4. Axes: distance in 


parsecs. 
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Model Bq\ spectrum, observed at time, t=4. 
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Model A 3 : X-ray spectra, observed at 6 = 45°, as a function of time 





